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We have introduced a new perturbative approach for t— J — V model where Hubbard operators are 
treated as fundamental objects. Using our vertices and propagators we have developed a controllable 
large- N expansion to calculate different correlation functions. We have investigated charge density- 
density response and the phase diagram of the model. The charge correlations functions are not 
very sensitive to the value of J and they show collective peaks (or zero sound) which are more 
pronounced when they are well separated (in energy) from the particle-hole continuum. For a given 
J a Fermi liquid state is found to be stable for doping 5 larger than a critical doping 5 C . 5 C decreases 
with decreasing J. For the physical region of the parameters and, for S < 8 C , the system enters in 
an incommensurate flux or DDW phase. The inclusion of the nearest-neighbors Coulomb repulsion 
V leads to a CDW phase when V is larger than a critical value V c . The dependence of V c with 5 
and J is shown. We have compared the results with other ones in the literature. 
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I. INTRODUCTION 

In the last years a large part of the solid state com- 
munity has been devoted to understand the physics of 
strongly electronic correlated models. The importance of 
this study is supported by clear experimental facts. High- 
T c superconductors^ and organic superconductors^ are 
considered, for instance, electronic correlated systems. 
They show that electronic properties are very different 
from those expected in usual metals. 

For better understanding of the physics of these sys- 
tems, it is very important to develop different methods 
for studying models for correlated electrons. The two di- 
mensional (2D) t — J model is probably one of the most 
simple models we should understand in the first place. 

The t — J model is the strong coupling version of the 
Hubbard modeli In spite of its simple form, there is 
no exact solution for this model until now and, many 
analytical and numerical 4 methods have been developed 
and the obtained results confronted to each other. 

In the t — J model double occupancy is forbidden 
and, the model can be written in terms of the corre- 
lated or projected operators, also known as Hubbard X- 
operators 5 . Hubbard operators verify complicated com- 
mutation rules which are very different from the familiar 
commutation rules for usual fermions and bosons. One of 
the common methods, used for avoiding this problem, in- 
troduces slave-particles^ in order to decouple the original 
A-operator. These slave-particles verify usual commuta- 
tion rules but, they are fictitious and sometimes it is not 
clear if the obtained results are genuine or artifacts from 
merely decoupling. On the other hand, the decoupling 
scheme introduces a gauge degree of freedom which re- 
quires a gauge fixing. The gauge fixing is a long discussed 
fundamental problem (see for example Ref.Q). 

Another problem in the treatment of the t — J model is 
the absence of any small coupling parameter suitable for 



a perturbative expansion. To deal with this problem, 
strong-coupling techniques have been developed. Be- 
tween them, what will be important for the present pa- 
per is the large- A expansion where N is the number of 
electronic degrees of freedom per site (see below). The 
large- A expansion has been used extensively in the con- 
text of slave-boson approach 7 (SBA) and, more recently, 
using Bayn-Kadanoff functional theory (BKF) in terms 
of A-operators^. 

On the basis of our Feynman path integral represen- 
tation for the t — J modeliS, we developed a large- 
N approach!! for the J = case ([/-infinite Hubbard 
model) . This method has the advantage of working with 
Hubbard operators as fundamental objects without any 
decoupling procedure. Then, we neither take care of 
gauge fluctuations nor Bose condensation as in the SBA. 

The method developed in Ref.^lJ was recently applied 
to describe electronic properties of quarter-filled organic 
molecular crystals^. For instance, in Ref.^2] we have 
calculated electronic self-energies S(fc, u>), spectral func- 
tions A(k,us) and the electronic density of states N(u) 
which are nontrivial calculations considering that they 
involve fluctuations, in a controllable approach, above 
the mean field. Our A(k, u>) and N(oj) were carefully 
compared^ with similar results obtained by Lanczos di- 
agonalization. Good agreement has been found for the 
behavior of the dynamical properties at J = 0.0. 

Our previous studies for J = 0.0 motivate us to extend 
our approach to the case of finite J. On the other hand, 
many physical problems play with a finite value of J and 
therefore, it is necessary to extend our approach to this 
case. This extension is the main topic of the present 
paper. 

As our method is new and we claim it can be used 
to simplify several calculations, we must show: a) how 
our method explicitly works and, b) comparison of the 
results with others in the literature in order to show the 
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confidence of the method. To satisfy a) and b) is one of 
the main purpose of the present paper. 

In Sec. II we develop the perturbative expansion and we 
give the new Feynman rules which are used in the explicit 
calculation in Sees. Ill and IV. In Sec. Ill we show results 
for the charge-charge correlation functions. In Sec. IV we 
study different instabilities of the model. In Sec. V we 
give the conclusions. 

II. PERTURBATIVE APPROXIMATION: A 
LARGE- iV APPROACH 

In this section, we present the large- ./V expansion for 
the t — J — V model in the framework of the path integral 
representation for Hubbard operators. We will extend 
our formalism of Ref.^J to the case of finite J. 

As in Ref.^3|, our starting point for developing the 
large- ./V approach is the path integral partition function 
written in the Euclidean form (it — ► r) , 

z = J vxf six™ + ^ xr i] s[xr' xr ° x x !° } 

x (sdetM AB )f exp (- J dr L E (X, X)) . (1) 

In Eq.(l), the five Hubbard V-operators X aa ' and V 00 
are boson-like and the four Hubbard V-operators V CT ° 
and X 0<T are fermion-like*. The spin index a is a = ± 
(up and down state, respectively). 

There is also a superdeterminant (sdetMAs) 1 = 
1/ (_j^oo)2 formed with the set of all the second class con- 
strains of the theory (see Refs.[n| and [13). 

The Euclidean Lagrangian L E (X,X) in (1) is: 

L E (X, X)= 1 -J2 ^ X *°+J r ^ + H(X) .(2) 

i,cr 1 

i 



On the basis of Hubbard V-operators, the t — J — V 
Hamiltonian is of the form: 



H(X) = Y, (tijXfX^ + h.c.) 

<ij>,<r 

+ \ E hi(xfxf-xrxf) 

<ij>;cr 

+ £ Viixrxf*' xr. (3) 

<ij>;(Ta' i.er 

Now, similarly to Ref.^j] we make the following 
changes in the path integral (1): 

a) We integrate over the boson variables X aa using 
the second S -function in (1). 

b) The spin index a = ±, is extended to a new index 
p running from 1 to N . In order to get a finite the- 
ory in the N- infinite limit, we re-scale the hopping tij, 
the exchange parameters Jy and the nearest-neighbors 
Coulomb repulsion Vij to ty/N, J%j/N and Vij/N, re- 
spectively. 

c) The completeness condition (X?° + J2 P Xf = N/2) 
can be exponentiated, as usual, by using the Lagrangian 
multipliers A^. 

d) The charge- like terms (the 3rd and 4th term of (3)) 
of the Hamiltonian, after extended to large- TV, are writ- 
ten in terms of X 00 using the completeness condition. 

e) We write V 00 and A in terms of static mean-field 
values and dynamic fluctuations; Xf° — Nr$(l + SRi), 
Aj = Ao + SXi. 

f) Finally, we make the following change of variables; 
f + — i y p ° f — 1 y° p 

By following the steps a-f, we find the effective La- 
grangian 



L eff ~ + 



/ tp 



(1 + SRi) 



<ij>,p 



(1 + SRi) 



Jij fipfip' fjp'fip 

2 ^ ~N (1 + SRi) (1 + SR,) 



Nrl J2 (Vij ~ Ui j )SR i SR j 



<y> 



(1 + SRi) 



SXi 



J ghost ■ 



(4) 



(sdetMAs) 2 leads to L g ] lost (Z) = Now, we treat the exchange terms Jy. These can be 

_ y\ Z \ ( l.) Zip when is written in terms of decoupled in terms of the bond variable A - through a 

p p \ + tt, j Hubbard- Stratonovich transformation, where A y - is the 
complex boson ghost field 



field associated with the quantity V &h . Finally, the Lagrangian (4) results 



Leff — ^ iyfipfip + fipfipj Q 

i,p 



fipfjp 



r 



Ay + /i.e. 



(5) 



At this point it is necessary to discuss similarities and 
differences between the SBA and our approach. 

a) From step f), we see that the fermions f p are pro- 
portional to the constrained X 0p operators. They are not 
associated with the spinons as in the SBA. 

b) From step e) , the field 5R is proportional to the real 
AT 00 operator representing the number of holes (empty 
sites). This is not associated with the holons as in the 
SBA. 

c) The bond variable Ay looks close to the valence 
bond variable of the SBA. However, A t j besides to be a 
function of the correlated fermions f p , it is also a function 
of SR through the denominator. 

Note that £ e // (Eq.(5)) contains several nonpolyno- 
mial terms. These apparent complications are the price 
we have to pay for working in terms of A-operators. Since 



the t — J — V model Hamiltonian is quadratic in the A's, 
the strong electronic interactions are contained in the 
commutation rules and the constraints. In the path inte- 
gral formulation, the information contained in the Hub- 
bard algebra was transfered to the effective theory. 

In the next sections we will show that present theory 
can go beyond a formal level and, the obtained results 
can be compared with others in the literature. 

Now, we write the Ay fields in term of static mean- 
field values and dynamics fluctuations A^ = A(l + r( + 
iAj), where rj = x,y and r( and correspond to the 
amplitude and the phase fluctuations of the bond variable 
respectively. 

To implement the l/N expansion, the nonpolynomial 
L e ff should be developed, as in Ref. 111. in powers of SR. 
Up to order 1/A the following Lagrangian is sufficient 



L eff = - 2 E (4/5 + fihp) (1 - siu + 5Rt) + ]r (t ijr j+f jp + h.c.) - M ^ ftpUi 1 - SR * + 6R i) 

i,p <ij>,p i,p 

OAT 

+ N r Q £ SX t SR, + ftpfiv^ ~ SIU) SX, + — A 2 £ [(rf ) 2 + (A'/) 2 ] 

i i,p it] 

- A E (f i pj P ' + f^Jip)[l-\m + SR j ) + l5R i 5R j + ^(5R^ + SR 2 j )} 

<ij>;P,P' 

- A E (tfpfjp'+f3 P 'fip)W+iAV)[i- Im+sRj)} 

<ij>,P,P' 

+ Nr% ]T (Vy - 1 -J ij )5R i 8R j - £ zj p (l - 5R t + 5R*) Z ip , (6) 

<ij> ip 

I 



where we have changed /i to /1 — Ao and dropped constant 
and linear terms in the fields. 

Looking at the effective Lagrangian (6), the Feynman 
rules can be obtained as usual. The bilinear parts give 
rise to the propagators and the remaining pieces are rep- 
resented by vertices. Besides, we assume the equation (6) 



written in the momentum space once the Fourier trans- 
formation was performed. 

In leading order of l/N, we associate with the N- 
component fermion field f p , connecting two generic com- 
ponents p and p' , the propagator 
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Feynman Rules 

Propagators 




P' 



FIG. 1: Summary of the Feynman rules. Solid line repre- 
sents the propagator G( ) (Eq.(7)) for the correlated fermion 
f p . Dashed line represents the 6x6 boson propagator 73( ) 
(Eq.(8)) for the 6-component field SX a . Note that the com- 
ponent (1, 1) of this propagator is directly associated with 
the X 00 charge operator. Doted line is the propagator B 
for the boson ghost field Z v . A. v a p (Eq.(9)) and repre- 
sent the interaction between two fermions f p and one and two 

bosons 5X a respectively. T„ p and F^ represent the inter- 
action between two ghost fields Z p and one and two bosons 
5X a respectively. 



8 i 

w(^») = -^r^ (7) 



which is 0(1). 

In (7), Ek = —2(tr + A)(cosk x + cosk y ), is the elec- 
tronic dispersion in leading order, where t is the hopping 
between nearest neighbors sites on the square lattice. 

The mean field values ro and A must be determined 
minimizing the leading order theory. From the complete- 
ness condition ro is equal to 5/2 where <5 is the hole doping 
away from half-filling. On the other hand, minimizing re- 
spect to A we obtain A = ^ ^ fc cos(k x )riF(Ek — /i) 
where uf is the Fermi function and N s is the number of 
sites in the Brillouin zone (BZ). 

For a given doping r5, the chemical potential /i and 
A must be determined self-consistently from (1 — 5) = 

w:T,k n F{Ek- n). 

We associate with the six component SX a = 
(5R , 5X, r x , r v , A x , A v ) boson field, the inverse of the 
propagator, connecting two generic components a and b, 



/ (4V -2J)rl(cos{q x ) 










cos(q y j) r 































7A 2 








^A 2 



(8) 



The bare boson propagator D^ ab (the inverse of 
Eq.(8)) is 0(1/N). As we will see D is renormalized 
to D by an infinite series of diagrams of 0(1/N). 



We associate with the N-componcnt ghost field Z p , the 
propagator, connecting two generic components p and p', 



VP 



-Sppi which is O(l). 



The expressions for three-leg and four-leg vertices are: 
a) 



(-1) ^K + ^'„)+M + 2A^ co8(k„-%)coa%; 1; -2 A cos{k x - -y) ; 



2 A co,[k h - %); 2 A sin(k x - %) ; 2A sin(k y - ^)) S pp ' 



(9) 



v 2" v "* 2 

represents the interaction between two fermions and one boson. 

b) A pp , which represents the interaction between two fermions and two bosons, is a 6 x 6 matrix where the only 
elements different from zero are: 



Krsh = (5K + + H + A E cos ( k v [cos% cos q j + cos'-^l}) 5™' (10) 



Irreducible Boson Self-Energy 



Af R rv = - A cos(k v - Q -^l) cos Q j 5* (12 

Af R A „=A Sl n(k v q -^l) cos q j 5* (13 

c) TP/ = (-1)(V . °> 0) °> 0) °) represents the 
interaction between two ghosts and one boson. 

d) is a 6 x 6 matrix, where Ijj^jj = 5 PP < and, th< 
others components are zero. It represents the interaction 
between two bosons and two ghosts. 

Each vertex conserves the momentum and energy and 
they are 0{1). In addition, in each diagram there is a 
minus sign for each fermion loop and a topological factor. 

Fig.l summarizes the Feynman rules. 

After identifying the propagators and vertices and the 
respective order of them, in the next sections, we will 
calculate different physical quantities. 

Before finishing this section, one remark is necessary. 
From the A-extended completeness condition we can see 
that the charge operator A 00 is 0(A), while the opera- 
tors X pp are 0(1). This fact will have the physical conse- 
quence that the 1/N approach weakens the effective spin 
interactions compared to the one related to the charge 
degrees of freedom. This is discussed in the next section. 

III. CHARGE CORRELATIONS 

In this section density-density correlations functions 
are calculated. It will be showed that the developed for- 
malism can be used in explicit calculation of different 
correlation functions and we will also compare our re- 
sults with others in the literature. 

The density-density correlation function is defined 

--11.13 



A; = ^£<7VAfAj V > (W) 
pp' 

Using J2 P X f P = N / 2 - X i° we nn d for D in the 
Fourier space 



Z?(q, LUn) = -AT(-) DsrsbH Wn) (15) 

The charge correlation, in 0(1), needs the calculation 
of all 0(1/A) contributions to DsR8R(q.,0J n ). From the 
Dyson equation, (-Dab) -1 = (£>(o)afc) _1 — II a 6, the dressed 
components D a b of the boson propagator can be found 



n = n (1) + n (2 > + n ,3) + n (4) 

~® -<y- * -Q- + -<>- + ----- 

FIG. 2: The four different contributions Ilj$ (i = 1, 2, 3, 4) to 
the irreducible boson self-energy H a b- 



after the evaluation of the 6x6 boson self-energy ma- 
trix Tl a b- Using the Feynman rules we may evaluate H a b 
through the diagramas of Fig. 2. Note that Dsrsr is the 
element (1,1) of the 6x6 dressed propagator. 

We note here an important difference with respect to 
the SBA. In SBA, only in leading order, the charge- 
charge correlation can be associated with the holon prop- 
agator. Beyond leading order, the convolution of holon 
propagators (which means the reconstruction of the X- 
operator A 00 ) is necessary. Meanwhile, in our case, it 
is not used any decoupling scheme and then, the differ- 
ent correlation functions are directly associated with our 
field variables. 

The density-density spectral function, which in princi- 
ple can be measure with electron energy loss scattering 1 *, 
is defined as D(q, uS) = —Im[D(q,ui)]. The imaginary 
part is taken as usual after performing the analytical con- 
tinuation iuj n = u> + ir\. 

In Fig. 3 we show the density-density correlation func- 
tion for the physical value J = 0.3 and for doping 
5 = 0.20. We plotted the densities for different q vec- 
tors in the BZ as a function of u. The q vectors and the 
physical parameters are the same to those used in the 
calculation of the densities in Ref . |15j . We used 77 = 0.1 
in the analytical continuation. 

Comparing Fig. 3 with Fig. 3 of Ref.flj|, we find a re- 
markable agreement between both methods in spite of 
the fact that the two approaches are very different. 

The density-density correlation function is nearly in- 
dependent of J. For example, comparing Fig. 3 (for 
J=0.3) with Fig.2 in Ref. pi (for J = 0) for momentum 
q = (71", 7r), we find a well pronounced collecti ve p eak (or 
zero sound) at w ~ 3.5. (Note that in Ref.|ll| lj is in 
units of t meanwhile in Fig. 3 it is in units of It in or- 
der to absorb the factor 1/N due to the re-scaling of the 
hopping term.) 

The fact that the density correlations are nearly in- 
dependent of J means that D(q, iu> n ) is dominated by 
charge fluctuations. As we mentioned above, the present 
formalism privileges charge over spin fluctuations. An- 
other consequence of this result is the absence, in O(l), of 
collective excitations (like magnons) in the spin suscepti- 
bility. The spin-spin correlation function is the electronic 
bubble with renormalized band due to correlationsii*^. 
Meanwhile there are collective effects in the charge sector 
in O(l), that appear in 0(1/N) in the spin sector. 
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FIG. 3: Density-density correlation functions for J = 0.3, 
V = 0.0 and 5 = 0.20 for different q's in the BZ. The density- 
density correlation functions contain collective peaks due to 
infinite diagram summation seen in Fig. 2. The collective 
peaks are well pronounced when they are above the particle- 
hole continuum (see for example q = (ir,0)). They are less 
pronounced when they are superimposed to the particle-hole 
continuum (see for example q = (2vt/5, vr/5)). U) is in units of 
2t and t is considered to be 1. 



FIG. 4: Phase diagram in the 8 C — J plane for V = 0.0. The 
solid line defines the stability border for the homogeneous 
Fermi liquid (HFL). For a given J the HFL is stable for 8 > 8 C . 
The thick dashed line around J ~ 0.5 separates the flux phase 
(FP) from the bond order phase (BOP). 



IV. INSTABILITIES 



The collective peaks are more pronounced when they 
are present well above the particle-hole continuum as for 
q = (tt, 7t). For momentum q — (2tt/5, 7r/5) the collective 
peak is superimposed to the particle-hole continuum and 
it appears broader. The broadening of the collective peak 
q = (tt, tt) is not intrinsic and it is only due to the finite 
value of r) we used in the analytical continuation. 

At this point we can compare our results with those ob- 
tained using exact diagonalizationi^. Despite the charge 
correlation function in Rcf. 16} was for doping S = 0.25 
and J = 0.4 and our calculation is for S — 0.20 and 
J = 0.3, there are some similarities and differences to re- 
mark. For example, Fig. 2 in Ref.^ij shows also a peak 
at q = (tt,it). However, a) this peak appears at larger 
energies than our collective peak, b) the peak at (tt, it) ob- 
tained in Ref.|l(j contains more structure than the ours 
and c) we must note also a difference at q = (tt, 0). In our 
calculation we found a collective-like peak at uj ~ 2 and 
a broad and small continuum at low energy. Lanczos 
diagonalization presents also this picture but, in addi- 
tion there is a peak at lower energy. As it was already 
pointed out by Khaliullin and Horsch in Ref. ( see a ls° 
Ref . [Isj ) . the inclusion of fluctuations beyond the mean 
field level is probably responsible for the differences listed 
in a), b) and, c). We think that the inclusion of fluctu- 
ations is more important at lower than at larger doping. 
In Ref. 13 j for doping 5 = 0.5, we found a better agree- 
ment between the peak position for q = (tt, tt) obtained 
by our approach and that obtained by Lanczos. 



A. Flux and Bond-Order phases 

The theory developed in previous sections defines a ho- 
mogeneous Fermi liquid (HFL) phase. In 0(1) the mean 
field solution (r , A) was independent of sites. Therefore, 
at A-infmite, we have free fermions with renormalizcd 
band due to correlations. 

In this section we will study the stability conditions for 
the HFL. 

In leading order, there are collective effects in the 
charge sector whereas there are not in the spin sector. 
Therefore, in principle, we expect instabilities in the 
charge channel. 

The HFL system is unstable when the static charge 
susceptibility Re(D(q,iui n — 0) diverges. The q vector, 
where the instability occurs (q c ), being the modulation 
of the new phase. 

Fig. 4 shows the phase diagram of the model for V = 0. 
For a given J, below a critical doping S c where the static 
charge susceptibility diverges at a given vector q in the 
BZ, the HFL is not the stable phase. 

In the limit J — > S c — > the HFL is stable for the 
whole doping range except at half filling where the system 
is an insulator because the band effective-mass tends to 
infinite. 

The instability is placed, for all J, on the border of the 
BZ. That is, q c = (1, x)tt or q c = (x, 1)tt where the pa- 
rameter x measures the degree of the incommensuration 
of the instability. 




0.1 0.15 

doping 5 C 



FIG. 5: Inconmensuration x vs S c for V = 0.0. x is defined 
as q c = (1,x)tt. q c is the q-vector where the instability takes 
place. The instability is incommensurate for all J except for 
J — > where x — > (q c — > (tt, tt)). At the separating border 
of the FP from the BOP there is a jump in the value of q c . 
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FIG. 6: Static charge susceptibility Re(D((\,iuj n = 0)) vs 
q = (tt, q) for J = 0.3, V = 0.0 and, for 5 = 0.15 and 5 = 0.14. 
For 5 = 0.14 (solid line) the curve looks similar to the case for 
S = 0.15 (dashed line). For S = 0.14 the HFL is not stable and 
there is no strong sign of the instability in the static charge 
susceptibility. In the inset we plot the det(D~^) (see text for 
discussions) vs q = (it, q). For 8 — 0.14, the det(D~^) clearly 
changes sign in a small region around q c ~ (it, 2.7). 



In Fig. 5 we plot the incommensuration 1 as a function 
of S c . For S c — ► 0, x — ► 1. Therefore, in the limit of 
J — ► 0, wc found a phase with commensurate order q c = 

(7T,7r). 

It is easy to overlook the instability just by looking at 
the static charge susceptibility. In Fig. 6 we plot the static 
charge susceptibility, along q = (n,q), for 5 = 0.15 and 
5 = 0.14 and, for the physical value J = 0.3. From Fig. 4 
we know that the HFL is not stable for doping 5 = 0.14 
whereas for 6 = 0.15 it is stable. Both curves in Fig. 6 
look similar and, for 5 — 0.14 there is no indication of the 
instability for q ~ q c because it occurs in a very narrow 
region of momentum near q c . In practice, this means 
that the set of q-points in Fig. 6 is not dense enough to 
localize the divergence of Re(D((\, iu> n = 0)). As we will 
see below, this is related with the fact that the instability 
is weakly coupled with the charge sector. 

To calculate -D(q, u> n ) we have to evaluate the inverse 
of (Dab)^ 1 . Therefore, for a better determination of the 
instability, we may look for the zeros of the determinant 
of . In the inset of Fig. 6 we plot the det(D~^) as a 
function of q = (tt, q). For 5 = 0.14 (close to the onset of 
the instability) the determinant changes sign in a narrow 
region around q c ~ (tt, 2.7) making the HFL unstable. 
Meanwhile, for 5 = 0.15 the determinant is positive for 
all q and the system is stable. A similar plot for 5 < 0.14 
shows, of course, a larger negative region for det(D~^). 

In order to characterize the nature of the new phase, 
we have studied the eigenvector corresponding to the zero 
eigenvalue of (Z? a b) _1 . 

For doping 8 < S c ~ 0.5 the eigenvector corresponding 
to the zero mode is ~ (0, 0, 0, 0, —1,1). Then, in the new 
phase, the 5th and 6th components of SX a (A x and A y ) 
are frozen. A x and A v are associated with the phase of 
the field A v and they take opposite values. Therefore, 



as the instability occurs near the momentum (7r,7r), this 
new phase is the well known flux phase (FP) 19,20 which 
opens a gap, with d-wave symmetry in the normal state. 

Recently, Chakravarty, Laughlin, Morr and Nayak con- 
sidered that the FP or DDW as a candidate to explain 
the physics of the pseudogap in underdoped cuprates2i. 

The first component of the eigenvector corresponding 
to zero eigenvalue is not zero exactly. It has a small value 
which means that the FP is weakly coupled to the charge 
sector. It means that the DDW is not completely hidden 
when it is incommensurate. Although Bragg peaks are 
predicted, they will show low intensity which would make 
their observation very difficult. This is the reason by 
which, any charge-probe is not very sensitive to show the 
FP instability. When q c is exactly (tt, tt) the eigenvector 
does not have any mixing with the charge sector and, the 
flux phase is fully hidden. 

For doping S > 5 C ~ 0.5 the unstable eigenvector is of 
the form ~ (0, 0, 0, 1, 0, 0). This state is known as bond- 
order phase (BOP)22iS. 

In the range of the studied parameters we did not find 
indications for phase separation. 

In Fig. 5, for J ~ 0.5, near the crossover from the FP 
to BOP, there is a jump in the incommensuration x. 

The agreement with other methods^S* 2 ^ is again re- 
markable in spite of the fact that the present approach 
is, a priori, very different from those ones. Taking into 
account that we constantly find similar or even the same 
results by means of different methods turn the results re- 
liable. For example, for the physical value J = 0.3, the 
HFL is unstable against a FP. This result is very robust 
and, many different methods agree that the stable phase, 
for small doping, is a FP for the physical region of the 
parametero 19 ! 20 ! 22 ! 23 ! 24 . 
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B. Charge Density Wave phase 



The inclusion of a nearest-neighbors Coulomb repul- 
sion V favors a charge density wave (CDW) state. In 
order to investigate this instability, we have included a 
finite value of V in the calculation. 

The study of a nearest-neighbors Coulomb repulsion 
in correlated models is important. In cuprates, there 
are indications^ showing the existence of V with a value 
close to J. If J favors superconductivity and, the value of 

V is large, it is reasonable to think that superconductiv- 
ity will be diminished due the presence of V. There are 
some analytical and numerical works studying the com- 
petition between V and J on superconductivity. Mean- 
while some papers indicate that superconductivity in the 
t — J model vanishes for V ~ J 9 , others indicate that 
superconductivity survives up to values of V >> 

On the other hand, the presence of V seems to 
be important for understanding the physics of organic 
materials^. Some peculiarities in the optical conductiv- 
ity (and also superconductivity) occur near the charge 
order and this picture can be interpreted by the com- 
petition between V and the kinetic energy in correlated 
models 2 ^. 

Looking at the Eq.(8) we see that V is only present 
in the element Dzi(l, 1). The Coulomb term enters in 

our approach multiplied by {5/2) 2 which means that the 
effect of V, at low doping, is strongly screened by the 
correlations. 

We note that in the SBA V can enter in different man- 
ners depending on the way chosen for the decoupling. 

In Fig. 7 we show the phase diagram in the V c — S plane 
for J=0 and J = 0.3. The curves show, as a function of 
6, the critical Coulomb repulsion V c where the CDW in- 
stability takes place. For V < V c the HFL is stable. This 
result is in agreement with numerical^ and analytical 
methods 29 . In Ref.|29j the authors use coherent poten- 
tial approximation and their results are close to the ours. 

A similar result was also recently found, at J = 0, 
in the triangular lattice in the context of the new low 
dimensional superconductor Na x CoO^&. 

For J = 0.3 and doping S = 0.20, the HFL is stable at 

V = 0. When V = V C ~ 1-25 the system enters in a CDW 
state. This V c is ~ 40% larger than for ,7 = and the 
same doping. There are two sources for this tendency. 

a) When J is finite, there is a contribution A to the 
effective hopping in Then, the system win kinetic 
energy when J is finite and, a larger V is necessary in 
order to localize the charges. 

b) As we see in the element (1, 1) of D7 Q y the effect of 

V is diminished when J is finite. The term 2 J in the ele- 
ment (1, 1) comes from the charge-like term JijXf a X° a 
of the pure t — J model, which is of the same form of the 
Coulomb term (see Eq.3). 

In order to identify the main source for the increas- 
ing of V c we looked for the CDW instability without the 
charge-like term in (3). Under this condition, for J = 0.3 
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FIG. 7: Charge density wave (CDW) stability line V c vs 6 for 
J = 0.0 (dashed line) and J = 0.3 (solid line). For a given S, 
the system is in a CDW state for V > V c . The value of V c 
increases with increasing J. For J = 0.3 the first presented 
doping is 5 — 0.15. (For J — 0.3 the HFL is not stable for 
S < 0.14 at V = 0.0). 

and 5 = 0.20, V c is V c ~ 1.0. This value is nearly the 
same than V c for J = 0. Therefore, b) is the main source 
for the increasing of V c with increasing J. 

In contrast to the FP and BOP, the CDW instability 
can be detected clearly by looking at the static charge 
susceptibility i?e(x(q, iuj n = 0)). 

In Fig. 8 we show, for J = 0.3 and S = 0.25, 
i?e(x(q, iu n = 0)) vs q = (q,q) for V = 1. (a little 
smaller than the critical value 1.1) and V = 0.5. In 
the figure, we see that approaching V c , the static suscep- 
tibility tends to diverge near q = (w, tt). For V = 0.5 
(V << V c ), the static charge susceptibility is flat. 

In this case the eigenvector corresponding to the zero 
eigenvalue is ~ (1,0,0,0,0,0). Then, the instability is 
mainly in the pure charge sector. 

From Fig. 8, in agreement with other methods^, the 
instability q- vector is always at q = (7r,7r). Therefore, 
in the CDW phase, the charges order themselves form- 
ing a checkerboard pattern. For doping different to the 
commensurate one 5 = 0.5 we must interpret the CDW 
instability as the system is charge-ordered but remains 
metallic. These results were obtained by different meth- 
ods and the nature of the instability for doping different 
to the commensurate one is still open. For the commen- 
surate doping S = 0.5 in Ref.[ll| we showed that the 
system is charge-ordered at V c ~ 0.6 however, we also 
find an insulator because the quasiparticle weight goes 
to zero exactly at this V c value. 

In order to compare the above results to the predictions 
of the BKF method, we calculated the critical V c using 
the formulation of Refs.0 @ for J = 0.3. For S = 0.20 
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FIG. 8: Static charge susceptibility for J = 0.3 and 8 = 0.25 
vsq= (q, q) for two different values of V, V = 1. and V = 0.5. 
V = 1. is near (from below) the value of V c = 1.1. In contrast 
to Fig.6, the CDW instability has a clear signal in the static 
charge susceptibility. For V = 1. there is a clear indication 
of a divergence at (n, n). For V = 0.5 (V << V c ), the static 
charge susceptibility is flat. 

and 6 = 0.60 we found V c ~ 1.25 and V c ~ 0.9 respec- 
tively. The agreement with our calculation is again very 
good in spite of the fact that the Coulomb term enters 
BKF in a different way. 



V. CONCLUSIONS 

From the path integral representation for the Hubbard 
X-operators we have developed a new perturbative the- 
ory for the t — J — V model. Our formulation is free of 
slave-particles and the X-operators are treated as funda- 
mental objects. 

We have extended the large- N formalism, recently de- 
veloped in Ref. to the case of finite J. 

To work directly with the X-operators has several ad- 
vantages upon the slave-particles methods. For example, 
in the slave-boson theory, the decoupling procedure intro- 
duces a gauge degree of freedom and then, an additional 
problem as the gauge fixing arises. This problem does not 
occur in our formulation. On the other hand, in the slave 
approaches, in contrast to our method, beyond a mean 
field level, a convolution between slave-particles propa- 
gators is necessary to reconstruct the original physical 
propagator. 

Fig.l summarizes the Feynman rules of the present ap- 
proach where the propagators and vertices are written in 
terms of the Hubbard operators X a0 , X 00 , etc. 

It is important to note that we can not read the inter- 
action vertices from the t and J terms of the Hamilto- 



nian. They arise from the nonpolynomial effective theory 
defined by L e // in Eq.(6). All the interactions caused 
by the Hubbard algebra (commutation rules and con- 
straints) were transfered to L e ff. 

The formalism developed in the present paper privi- 
leges charge over spin fluctuations. While there are col- 
lective effects in the charge sector in 0(1), those appear 
in 0(1/ N) in the spin sector. 

We have studied charge correlations functions and we 
have investigated the role of collective effects. This study 
shows the presence of collective peaks (or zero sound) 
which are more evident when they are separated from the 
particle-hole continuum. On the other hand, the collec- 
tive peaks become broader when they are superimposed 
to the particle-hole continuum. One important charac- 
teristic of the charge correlation function is the fact that 
this is not strongly dependent on the value of the ex- 
change interaction J. 

At large N, the theory can be described as a homoge- 
neous Fermi liquid with renormalized band due to cor- 
relations. We have studied the stability of this Fermi 
liquid phase. For a given value of J, the Fermi liquid 
phase is stable for doping 8 > S c . The critical doping 
S c decreases with decreasing J. For doping below S c the 
system enters in a flux or bond order phase depending 
if J < 0.5 or J > 0.5 respectively. These instabilities 
are weakly coupled to the charge sector which means 
that any charge-probe is not an efficient test to detect 
them. One important characteristic of these new phases 
is that they are incommensurate with a modulation vec- 
tor q = (l,x)ir where the incommensuration x tends to 
1 when J goes to zero. 

It is important to remark that in agreement with other 
theories, for low doping, the Fermi liquid is unstable 
against a flux or <i-density wave phase for the physical 
region of the parameters. 

We also have investigated the role of a nearest- 
neighbors Coulomb repulsion V on the stability condi- 
tions of the Fermi liquid. When V is larger than a critical 
value V c , the system enters, at a given doping, in charge 
density wave state. The value of V c increases with in- 
creasing J. We have identified in the charge-like term of 
the pure t — J model the main reason for this increase. 

We have continuously compared our results with sim- 
ilar ones in the literature. The agreement of our results 
with those obtained by other methods give confidence to 
the results and the approach of the present paper. 

At leading order, our formalism is in agreement with 
the slave-boson approach. However, at the next to lead- 
ing order (which is necessary to calculate dynamical 
properties) the differences between the two formulations 
are not yet completely established. We think that the 
complications that appear in the slave-boson method be- 
yond mean field level as gauge fixing, Bose condensation, 
regularizing factors^ are not good for the advance of the 
field. In our case we do not have these problems and 
we think that our approach can be useful to go beyond 
the mean field level. With the formulation for finite J in 
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hand, we expect to continue in this direction as we did The authors thank to C. Genz, L. Manuel, J. Merino 
for J = in Ref.^2- an d R- Zeyher for valuable discussions. 
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